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Abstract 



Duality relations are obtained for correlation functions of the g-state Potts 



C^ . model on any planar lattice or graph using a simple graphical analysis. For 
the two-point correlation we show that the correlation length is precisely the 

Vh ■ surface tension of the dual model, generalizing a result known to hold for the 

^' 

P^ . Ising model. For the three-point correlation an explicit expression is obtained 

o, 

, ^,. relating the correlation function to ratios of dual partition functions under 



fixed boundary conditions. 



I. INTRODUCTION 

It is now well-known |]l]^ that the correlation length of the Ising model in two dimen- 
sions is precisely the surface tension of the dual lattice. It is also known by folklore (see, 
for example, [^) that a similar duality exists for the Potts model Q]. However, detailed dis- 
cussion of the correlation duality for the Potts model has yet to appear in the literature. In 
view of the important role played by such duality relations in discussions of the equilibrium 
crystal shapes [§|,§[, a definitive understanding of this subject matter is clearly needed. 

Here we take up this question and consider more generally the duality for the n-point 
correlation function. On the basis of a simple graphical analysis, we derive duality relations 
for the two- and three-point correlation functions. Particularly, we establish that the corre- 
lation length in the large lattice limit is precisely the surface tension of the dual model, thus 
generalizing the known Ising result, and that the three-point correlation function is given 
by similarly defined expressions in the dual space. Our analysis, which does not concern 
with the internal structure of the underlying lattice, is quite general and applies to the Potts 
model with arbitrary edge-dependent interactions and/or on any planar lattice or graph. It 
can also be extended to higher correlations. 

Consider a g-state Potts model |^ on a two-dimensional lattice L with a free boundary. 
An example of L is the rectangular 4x6 lattice shown in Fig. 1(a), but more generally L can 
be any planar lattice or graph which does not need to be regular. The site i of L is occupied 
by a spin ai which can take on values 1,2, ...,g. Two spins at sites a and (3 and in states 
(Ja and (jp interact with a Boltzmann factor exp[i^o/3^(o"a, '^p)]- Let i,j,---,n denote any n 
sites on the boundary. We introduce the n-point correlation function as the probability that 
the spins at sites {i,j, ■ ■ ■ ,n} are in the respective states {a, a',- ■ ■ , cr*^"~^^}, 

P„(a, a', ■ ■ ■ , a("-i)) =< (5(a„ a)5ia,, a') ■ ■ ■ (5(a„, a^^-'^) >, (1) 

where < ■ > denotes statistical averages. It is sometimes convenient to consider the corre- 
lation function (see, for example [^^ for n = 2) 



Tn{ai,aj,-- ■,an) = < q"^ ^5{ai,aj)5{aj,ak) ■ ■ ■(5(cr„_i, (t„) - 1 > 

= g"P„(a,a,---,a)-l, (2) 

a quantity which vanishes identically if the n spins {(Xj, ctj, a^ • • ■ , cr„_i, a„} are completely 
uncorrelated. In the case of g = n = 2, for example, one writes 26{ai, ctj) = 1 + aiaj where 
aij = ±1, then (||) becomes the usual expression T2{cri,aj) =< aiaj > for the Ising model. 
It is clear that the correlation function P„ is more general than r„. 

Our analyses is based on the repeated use of a fundamental duality relation for the Potts 
model [0. Construct L*, the dual lattice (graph) of L, by placing spins in the faces of L 
including the exterior one. In Fig. 1(a), for example, the faces of the dual lattice are denoted 
by crosses and the letter s, so L* has A^* = 15 + 1 = 16 sites. Edges of L* (not shown in Fig. 
1(a)) bisect, and are in one-one correspondence with, edges of L. Let K^^ be the interaction 
of the edge bisecting the interaction K^is- Let Z{{Kaf3}) and Z*({_K'*^}) be the respective 
partition functions on L and L*. Then one has the following fundamental duality relation 

i 

Z{{K^f,})=qCZ%{K:^}) (3) 

where Kap and K^^ are related by 

{e^'^" - l){e^*-p - 1) = q, (4) 

andC = g-^*nedges(e^"''-l). 

II. THE TWO-POINT CORRELATION FUNCTION 

Let i and j denote any two sites on the boundary of L as shown in Fig. 1(a). Further let 

Zu = Z(ai = Gj = 1), Zi2 = Z(oi = 1, aj = 2) (5) 

be the partition functions with ai and aj fixed in definite states 1 and/or 2, Then, by 
symmetry we have 



Z{{K^p}) = qZu + q{q - l)^i2. (6) 

Clearly, we also have ^2(1; 1) = Zn/ Z and P2(l; 2) = ZyijZ. 

Let Z\^ denote the dual partition function with the spin in the exterior face in a definite 
state, say, s. Alternately, from Fig. 1(a), we see that Z\y can also be regarded as the dual 
partition function with all boundary spins interact with the exterior spin fixed at state s. 
Then we have 

Z\{K^f,})=qZl,. (7) 

Substituting (|]) and (0) into the duality relation (|), one obtains the relation 

Zn + (g - 1)^12 = qCZl^. (8) 

To obtain a second relation relating Zx\ and Z12, we apply the duality relation (§) to a 
lattice, or graph, constructed from L by adding one additional edge connecting sites i and 
j with an interaction K^ a situation shown in Fig. 1(c). Clearly, the respective partition 
functions for the new lattice and its dual are now 

Z = ge^Zn + g(g - 1)^12 (9) 

and 

Z* = qe^'Zl, + qiq~l)Z*, (10) 

where ZI2 is the partition function of the dual under the boundary condition such that the 
spins in the two exterior faces are in definite states s 7^ s' as shown in Fig. 1(b). Here, K 
and K* are related by (H). Substituting (P) and ([T0|) into (H) and noting that A^* is increased 
by 1 in the new lattice, one obtains a second relation 

e^Zn + (g - 1)^12 = C(e^ - l)[e^*Zi*i + (? ~ 1)^12]- (H) 

We now solve (^ and (|TT|) for Zn and Z12. This leads to, after using (^ to eliminate 



Zn = C[Zl, + {q- l)Zl,] 

Z^2 = C{Zl^-Zl^). (12) 



Note that the interaction K introduced to facihtate calculations does not enter (|T2|). The 
expression (|T^ now leads to the desired expression for the correlation function. Particularly, 
using the identity Z = q^CZ^^, we obtain 






^11 



(13) 



and 



r2(a.,a,) = (g-l)(||). (14) 

For the Potts model the surface tension r is defined by 

where D is the distance between sites i and j, the two points where the boundary spin 
state changes from s to s' in Fig. 1(b). Our result ( [I^ ) now relates the correlation function 



on L to the surface tension on L* . Particularly, the known exponential decay e~^l^ of the 
two-point correlation r2 above the transition temperature Tc [0, where ^ is the correlation 
length, leads to the identity ^^^ = r, where r is the surface tension of the dual model below 
Tc. This generalizes the corresponding q = 2 result of the Ising model [|l],0 to all q. 

III. THE THREE-POINT CORRELATION 

Consider any three sites i,j, fc on the boundary of L as shown in Fig. 2. Here, to 
emphasize the generality of our consideration and the fact that the internal structure of L 
does not enter the picture, the lattice is shown symbolically as a shaded region. Let 

Z„„>„n = Z{pi = a, Oj = a', aj = a"). a, a', a" = 1, 2, ..., q (16) 

denote the partition function with ai,aj,ak fixed in definite states. Then, in analogous to 
(^) the partition function can be written as 



Z = qZui + q{q - l)[^2ii + ^121 + ^112] + qiq - l)(g - 2)^123. (17) 

In a similar manner, let Z*^,^„ be the partition function of the dual model under a boundary 
condition such that all boundary spins between sites i and j of L interact with a fixed spin 
state s", all boundary spins between j and k interact with a spin state s, and all boundary 
spins between k and i interact with a spin state s', a situation shown in Fig. 2. Then, the 
dual partition function can be written as 

Z* = qZl,,. (18) 

Substituting (0) and (plSf) into (|), one obtains the relation 

Zni + (g - l)(^2ii + ^121 + ^112) + (g - i)(g - 2)^123 = qCZl^^. (19) 



Next we modify L by connecting sites j and k with a new edge and apply the result (|12|). 
In the first line of (|T2|) , one has, by definition, 

^11 = 2^ -^jii = ^111 + (g — 1)^211 

i=l 
^7* /7* 

^11 — ^111 

^12 = ^122 = ^211- (20) 

Substituting (^) into (|12|), one obtains a second relation 



Zni + (g - 1)^211 = C[Zl,, + {q- l)Z;,,]. (21) 

In a similar manner by connecting sites {k,i} or {i, j} with a new edge, we obtain, respec- 
tively, 

Zni + (g - 1)^121 = C[Z*ii + (g - l)Z*2i] (22) 

Zni + (g - 1)^112 = C[Z*,i + (g - l)Z*i2] (23) 

Finally, we apply the duality relation (^ to a lattice constructed by introducing 3 new edges 
connecting sites i,j or /c to a common new point with interactions K as shown in Fig. 2. 
This process increases N* by 2 and yields the duality relation 

6 



e^^ + g - l)Zni + {q - l)(e'^ + e^ + g - 2){Z2n + ^121 + ^ 



112i 



-l)(g-2)(3e^+g-3)Z 



123 



C(e^ - iyq-'[e'^ Z*,,, + {q - l)e^ (Z*n + ^1*21 + ^1*12) 
+ (g-l)(g-2)Z*23]. 



(24) 



From (|19]), (gT]), (||), (||), and (|2|) and using (|) to eliminate K*, we solve Zm, Z211, 
Z121, Z112, and Z123 in terms of the dual partition functions Z*^,^„. This leads to 



Zni = {C/q)[Z*,, + {q- l)(g - 2)Z*,, + {q - l)(Z*n + ^121 

^123 = (C/q) [^111 + 2^1*23 - ^211 - ^121 - ^112] 

^211 = iC/q)[Zl,, -{q- 2)Z*23 + {q - l)^2*ii " ^1*21 " ^1*12] 
Z121 = {C/q)[Zl,, - (g - 2)Z*23 - ^211 + (? " l)^i*2i " ^112] 
Z112 = {C/q)[Zl,, - (g - 2)Z*23 - ^211 " ^1*21 + (? - 1)^1*12]- 

This result is again independent of the parameter K used in the evaluation. 
Using (H) and the identity Z = CqZ* = Cq'^Z^^-^^, we finally obtain 



^1*12)] 



P3{cr,a',a' 



1 + 2pi23 - (P211 + P12I + P112) + q{P211 - Pl23)Sa' 



+ q{Pl21 - Pl23)Sa",a + g(Pll2 " Pl23)Sa,a' + q^Pl235u,a'5a' ,a" 



and 



(25) 



(26) 



T^{ai,aj,ak) = (g - 1) 



P112 + P121 + P211 + (g - 2)p 



123 



r2(o-i, (Tj) + T2{(jj, Ok) + r2(crfe, (Ti) + (g - l)(g - 2)pi23. 



(27) 



Here, 



Ps 



^ss's" I "^111 



(28) 



is the ratio of dual partition functions which can be interpreted as appropriately defined 
surface tensions, and we have used (g — l)pii2 = r2(o"j, aj), etc. Since Z* = qZ^^^, the ratio 



Pss's" gives also the probability that the particular boundary condition {s, s', s"} of Fig. 2 
occurs in the dual partition function Z*. It is readily verified that the identity 

Y^Ps{a,a',a") = P2{a,a') (29) 

a" 

is satisfied. 

IV. SUMMARY 

In summary, we have presented a graphical analysis of duality relations leading to explicit 
expressions for correlation functions of the g-state Potts model on any planar lattice or graph. 



For the two-point correlation our result (|T4D generalizes a relation previously known for the 



Ising model. For three-point correlations, our result (p6D is new and can be used to compute 
any three-point correlation. In all cases the correlation functions are found to be given by 
ratios of dual partition functions under fixed boundary conditions, which can in turn be 
interpreted as appropriately defined surface tensions. Our consideration can be extended 
in a straightforward fashion to higher correlations Pn and r„ for n > 3, and to the more 
general (NajNp) model Jll| including the A'^^ = A^/j = 2 Ashkin- Teller model |ll2 |. 
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FIGURES 
FIG. 1. (a) A 4 X 6 lattice L. The dual spins are denoted by crosses and the letter s. (b) The 
boundary condition {s / s') for Z^2- ('^) The addition of one edge to L connecting sites i and j. 

FIG. 2. The addition of 3 edges to L (the shaded region) connecting sites i, j, and A; to a 
common point. 
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